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Abstract 

A quantum relativistic model is proposed to represent the confining potential 
of cc and bb mesons by using the Klein Gordon oscillator. We have solved 
the two-body problem obtaining the eigenvalues and the corresponding mass 
spectra, which fit very well with experimental data. We compare the results 
of our model with others obtained from different quark binding potentials. 
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I. INTRODUCTION 


The study of bound sates of quark antiquark systems has been performed with different 
formalisms in the last two decades. These systems are a profitable testing ground for new 
ideas and methods to explain confinement and strong interaction dynamics. In particular, 
the mass spectrum for charmed mesons has largely been measured and modeled. Some 
time ago, Kang and Schnitzer have proposed a model to simulate confinement between 
the constituents of mesons. These authors considered a linear potential (LP) in the Klein 
Gordon equation. In this way they obtained the mass spectrum for charmonium. Later 
on. Ram and Hasala @] followed this procedure with a quadratic potential (QP). There 
are also some works in the literature, devoted to the study of quarkonium systems from a 
non-relativistic point of view, with different potentials. As an example of that Eichten et ah 
P] have studied the charmonium family with linear plus Coulomb interaction as a conhning 
potential. Furthermore, Quigg and Rosner [|] have studied the bound states of quarkonium 
with logarithmic, and power law potentials. 

In this work we introduce an alternative relativistic quantum mechanical model which 
includes the Klein Gordon Oscillator in order to obtain the mass spectra for bound states 
of bb and cc mesons. We present a scheme that considers quarks as spinless particles, in the 
same spirit of ref. [j^ and p, which reproduces with high precision the charmonium and bot- 
tonium mass spectra. The harmonic oscillator (HO) is a natural candidate to simulate the 
quark binding. It is one of the simplest potentials, well behaved, which takes into account 
asymptotic freedom and confinement. The Klein-Gordon oscillator (KGO) was introduced 
by Bruce and Minning p. This interaction gives the usual HO in the non-relativistic limit. 

In this frame we have solved the two-body problem. This alternative relativistic treat¬ 
ment is very suitable to study bound states due to its simplicity. We compare our results 
for charmonium, with those obtained for linear and quadratic potentials in the Klein Gor¬ 
don equation P p. Moreover, we have included in our analysis the values for the mass 
spectra for cc and bb mesons that we have obtained by considering the treatment for two 
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body systems in a quantum relativistic frame. This treatment was presented in a paper by 
Arshansky and Horwitz |p. There, it is derived the mass spectra and the wave functions for 
the 4-D harmonic oscillator (4D) as a quantum relativistic two-body bound state, but the 
authors have not £t their mass spectrum with the systems we are interested on, something 
we shall do below. 

The paper is organized as follows. In Section II we comment on the Klein-Gordon oscil¬ 
lator formalism. Section III is devoted to describe the model we propose. In Sect. IV we 
compare our results for charmonium with those obtained in refs. [|l| p|, and with those that 
we have obtained considering the formalism in ref. . We also present a table for the results 
of bottonium masses within KGO and 4D models. Finally, we sketch our conclusions. 

II. THE FORMALISM 

In 1989 Moshinsky and Szczepaniak proposed a new type of interaction in the Dirac 
equation. This potential is linear in both, momenta and coordinates, and the corresponding 
equation was named “Dirac Oscillator” because in the non-relativistic limit, the harmonic 
oscillator is obtained. Later on, this kind of interaction was also introduced in the Klein- 
Gordon (K-G) equation by making the minimal substitution: 

P —> F — irri'yfl.Q (1) 

where m is the mass of the particle and P = fjp , Q = fjq are respectively, coordinate and 
momentum operators. The matrix D is chosen to be 3x3 with coefficients Clij = The 

constants uji (with i = 1, 2, 3) are the oscillator frequencies along the x,y,z axes. The 
matrix 7 satishes the anticommutation relations: 

{7,h} = 0 f = 7]^ = l. 

Thus, the corresponding K-G equation (in natural units) for the wave function T [q, t) has 
the following form : 
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(jf + rn^q.Vt^.q + m^trVt + '^{q, t). 


In order to solve it, it is convenient to recast the K-G equation, and some ways to do this are 
shown in ref. |@]. One of these posibilities is the Sakata-Taketani approach |P] characterized 
by the equation: 




p(r3 + iT2)p 

2 m 


+ mrs } <h 


where r* are Pauli matrices and $ has the form 

* 

Uy 

This two component wave function is related with T as follows 


I + 

1 X = (^ - /2. 


( 2 ) 

(3) 

(4) 


In order to obtain the oscillator-like equation, in the isotropic case where a;i = a ;2 = c^s, the 
minimal coupling 


p —> p — imoJTir 


must be introduced in the K-G equation (H) giving: 

• ^ ^ f {v~ imu!Tif}(T3 + iT2){p — 1 ^ 


which yields 


□ -l-2ma;r — 3ma; +m T = 0 
The corresponding energy spectrum p] is 

E'^ = 2 muj {Ni + N 2 + iVa) -|- 


(5) 

( 6 ) 

(7) 

( 8 ) 


with iVi, N 2 , N 3 = 0,1... 

On the basis of the results presented in this summary, our model can be introduced. 
Our purpose is to work out the two-body problem in the frame of the formalism we have 
commented. 
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III. THE MODEL 


We propose that the above KG oscillator could model a quark antiquark conhning po¬ 
tential taking into account some further considerations. In the case of one body KGO, we 
note from (^) that the oscillator interaction has its origin at 0. If one desires to change the 
origin to rg, the replacement r — tq must be done. We are interested in studying the case of 
two interacting particles with the same mass. We call with subindex 1 and 2 the magnitudes 
corresponding to each particle. We propose a slight different ansatz when compared to (i- 
It reads, 


Pi 


Pi = Pi - irn^Txifi - rl) 


(9) 


P2^ fi =P2- im^Ti{r2 - f[), 


( 10 ) 


giving rise to two coupled equations: 


and 


• ^ as 
uti 


dt 2 


Pijrs + iT2)pi 
2 m 


P2{t 3 + iT2)p2 
2 m 


mT 3 


$ 


( 11 ) 


mra 


d). 


( 12 ) 


Summing the last two expressions, one obtains an equation which corresponds to the 
case of two particles: 

Pi in iT2)pi , pl'(r3 iT2)p2 


.9 -9 

uti ut2 


2 m 


+ 


2 m 


2 mr3 




(13) 


In quantum mechanics, the study of a two body system could be thought of as a free particle 
plus a particle in a central potential, which depends on the relative coordinate r = f\ — f 2 - 
In order to solve our system we have made the following change of variables: 


T = 


ti + t 2 
2 


t = ti-t2 


(14) 


and 
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R = 


n + r2 


r = n — r 2 - 


(16) 


The momenta pi and p 2 can be written as follows: 


^ P ^ 
Pi = l^+P 


(16) 


P2 = J - P- 


Shifting now to these new variables, eq. ( 0 ) gives: 




_|_ _|_ 4 , 7 ^ 2 ^ $(P,p, T). 


This equation can be solved by separation: 


(17) 


(18) 


‘i(R,f,T) = MR,T)Mr) 


(19) 


giving rise to the following two equations 

( 


dT^ 


+ P^ + E^\ 01 (P,T) = (□ + P')0i(P,T) 


and 


— Sfiuj + 4/r^ — ^ ) 02(^0 = 0 


( 20 ) 


( 21 ) 


where m = 2fi. We are particularly interested in solving the last equation, since the first 
one, has of trivial solution. In spherical coordinates this equation takes the form: 


Pr^ + ^ + - 3/ia; + 4/i^ - ^ ) 02(r) = 0. 




( 22 ) 


The first term can be written as: 


.( d V 

Pr = -M W + “ 

\or r 


(23) 


This Hermitian operator is the radial component of the linear momentum. There are two 
independent eigenfunctions. In analogy with non-relativistic quantum harmonic oscillator. 
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only one of them is physically acceptable because the other diverges at the origin. So, the 
eigenvalues and eigenfunctions are: 


- 16/1^ = 8/ra;(2n + /) n = 0,l... / = 0,1... (24) 

<h(r, 9, 0) (25) 


and An,I is the normalization constant: 


2(/ra;)'+ir(n + 1) 
7(12 + 1 + I) 


(26) 


Note that the states are degenerated for particular values of n and 1. The same situation 
appears in the 4D model. 

The solution we have previously derived enable us to compute the mass spectrum for 
quarkonium, as we will show below. 


IV. RESULTS AND CONCLUSIONS 

Let us now compute some specihc values for the mass spectra for charmonium and bot- 
tonium. We have considered a two free parameter set: rriq and oj. For the cc family rric = 
1.5 GeV and oj ~ 0.385 GeV. Instead, for bottonium, rrib = 5.17 Gev and oj ~ 0.12 GeV. 
The values for oj were obtained by adjusting the rms radius of the ground state P, 

Our results are shown in Tables I and II. Table I displays the values for the masses of 
charmonium mesons within different relativistic models, for particular quantum numbers n 
and 1. The explicit values for the 4D model arise from the computation of charmonium 
spectrum in such a general scheme. The other two known values, LP and QP, are taken 
from ref. With bold typos (KGO) are presented the values we have obtained with our 
model. From Table I, we can see that our scheme produces acceptable results, without any 
further assumption. 

In Table II, results for bb mesons are presented. Here, only the 4D model can be con¬ 
fronted with our one, because, as far as we know, the LP and QP treatments were not 
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applied to this mesons. To do such a comparison we choose the mass for bottom quark to 
be 5.16 GeV, which is bigger than the experimental value. This value for the mass was used 


in p] yielding lo ~ 0.113 GeV. Thus, we can expect an overestimation of the resulting mass 
spectrum, due to such a choice. Anyway, taken the b mass heavier than the experimental 
value, the masses for bb mesons predicted from these two models are very close to the ex¬ 
perimental data. We again get agreement between 4D and KGO. 

Summing up, we have presented a simple quantum relativistic scheme to take account 
the conhning potential of quarkonium. In this frame we have solved the two-body problem 
obtaining the corresponding eigenfuntions and eigenvalues. We have made a comparison 
between our model and other relativistic models. This family of mesons entails for a restric¬ 
tive arena in where comparison can be made, because plenty of models exist with agreement 
among them and experimental data. Thus, this can be thought of as a hrst test of any 
alternative setting. Our results are in good agreement with previous theoretical predictions 


and with experimental data (available in |^). This fact allows us to conclude that our 
treatment reproduces with a high degree of precision the masses of charmonium and botto- 
nium family. In addition, we stress that our model is enterely analytical. 

We expect that our simple model could be extended to other situations in the framework 
of quark models. 
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TABLES 

TABLE I. Particular values for cc mass spectrum are shown. In KGO and 4D models, we have 


considered rric = 1.5 GeV and u) = 0.385GeV. 



Relativistic Model 

1=0 

1=1 

1=2 

1=3 

l=A 


LP 

3.105 

3.456 

3.761 

4.034 

4.286 

n=0 

QP 

3.179 

3.946 

4.193 

4.436 

4.675 


4D 

3.531 

3.844 

4.133 

4.404 

4.659 


KGO 

3.000 

3.363 

3.691 

3.991 

4.271 


LP 

3.696 

3.964 

4.215 

4.450 

4.673 

n=l 

QP 

3.695 

3.946 

4.193 

4.436 

4.675 


4D 

4.133 

4.404 

4.659 

4.901 

5.131 


KGO 

3.691 

3.991 

4.271 

4.533 

4.781 


LP 

4.169 

4.397 

4.616 

4.826 


n=2 

QP 

4.186 

4.426 

4.662 

4.895 



4D 

4.659 

4.901 

5.131 

5.351 

5.563 


KGO 

4.271 

4.533 

4.781 

5.017 

5.242 


LP 

4.580 

4.783 

4.980 



n=3 

QP 

4.656 

4.886 

5.114 




4D 

5.131 

5.351 

5.563 

5.767 

5.964 


KGO 

4.781 

5.017 

5.242 

5.458 

5.666 
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TABLE II. Particular values for bb mass spectrum are shown. We considered rrii, = 5.17 GeV 


and uj = 0.113GeV. 



Relativistic Model 

1=0 

1=1 

1=2 

1=3 

1=4 

n=0 

4D 

10.358 

10.667 

10.795 

10.922 

11.047 


KGO 

10.340 

10.452 

10.564 

10.674 

10.782 

n=l 

4D 

10.795 

10.922 

11.047 

11.171 

11.293 


KGO 

10.564 

10.674 

10.782 

10.890 

10.997 

n=2 

4D 

10.047 

11.171 

11.293 

11.415 

11.534 


KGO 

10.782 

10.890 

10.997 

11.103 

11.208 

n=3 

4D 

11.293 

11.415 

11.534 

11.653 

11.770 


KGO 

10.997 

11.103 

11.208 

11.311 

11.414 
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